In [1]:
from nilearn import datasets, plotting
from netneurotools.datasets import fetch_schaefer2018
import numpy as np
import nibabel

/home/mikef/ENTER/lib/python3.8/site-packages/nilearn/datasets/__init__.py:86: FutureWarning: Fetchers from the nilearn.datasets module will be updated in version 0.9 to return python strings instead of bytes and Pandas dataframes instead of Numpy arrays.
  warn("Fetchers from the nilearn.datasets module will be "
In [2]:
## get schaefer:
fname_schaefer = fetch_schaefer2018("fsaverage5")['400Parcels7Networks']
schaefer_lh = nibabel.freesurfer.io.read_annot(fname_schaefer.lh)[0]
schaefer_rh = nibabel.freesurfer.io.read_annot(fname_schaefer.rh)[0]

## get mmp in fsaverage:
## from https://figshare.com/articles/dataset/HCP-MMP1_0_projected_on_fsaverage/3498446
mmp_lh_fs = nibabel.freesurfer.io.read_annot('in/lh.HCP-MMP1.annot')[0]
mmp_rh_fs = nibabel.freesurfer.io.read_annot('in/rh.HCP-MMP1.annot')[0]

## get mmp in fsaverage5:
mmp_lh_fs5 = nibabel.freesurfer.io.read_annot('in/lh.HCP-MMP1_fsaverage5.annot')[0]
mmp_rh_fs5 = nibabel.freesurfer.io.read_annot('in/rh.HCP-MMP1_fsaverage5.annot')[0]

## get mmp in fs_lr
mmp_lh_fslr = nibabel.load('in/HCP-MMP_RelatedParcellation210.CorticalAreas_dil_Colors.32k_fs_LR_L.label.gii')

## get underlays:
fsaverage = datasets.fetch_surf_fsaverage('fsaverage')
fsaverage5 = datasets.fetch_surf_fsaverage()
fs_lr_lh = nibabel.load('in/Q1-Q6_RelatedValidation210.L.very_inflated_MSMAll_2_d41_WRN_DeDrift.32k_fs_LR.surf.gii')
In [3]:
plotting.view_surf(
    fsaverage5['infl_left'], schaefer_lh, cmap = "prism", symmetric_cmap = False,
    title = "Schaefer 400 lh in fsaverage5"
)
Out[3]:
In [4]:
plotting.view_surf(
    fs_lr_lh.agg_data(), mmp_lh_fslr.agg_data(), cmap = "prism", symmetric_cmap = False,
    title = "MMP lh 400 in fs_lr"
)
Out[4]:
In [5]:
plotting.view_surf(
    fsaverage['infl_left'], mmp_lh_fs, cmap = "prism", symmetric_cmap = False,
    title = "MMP lh in fsaverage"
)
Out[5]:
In [6]:
plotting.view_surf(
    fsaverage5['infl_left'], mmp_lh_fs5, cmap = "prism", symmetric_cmap = False,
    title = "MMP lh in fsaverage5"
)
Out[6]: